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Abstract Instabilities in magnetic fields wound up by differential rotation 
as reviewed in Spruit (1999) are discussed with some detail and new develop- 
ments added. In stellar models which include magnetic torques, the differential 
rotation tends to accumulate in the gradients in composition. In view of this, 
instability in a /i-gradient is studied in more detail here, resulting in the detec- 
tion of a second instability. Its relevance for angular momentum transport is 
uncertain, however, since it requires high horizontal field gradients and would 
not operate near the pole. Finally, the possibility is discussed that magnetic in- 
stability in a /^-gradient will lead to layer formation: the gradients breaking up 
into small steps of uniform composition and rotation rate. This would enhance 
the angular momentum transport across inhomogeneous zones, and decrease 
the rotation rates of the end products of stellar evolution. A recent astro-ph 
submission by Dcnniscnkov and Pinsonncault proposing a modification of the 
instability conditions is shown to contain a mathematical error. 

1 Introduction 

1.1 Strong fields and weak fields 

Magnetic fields are potentially the main mechanism for transport of angular 
momentum in stars, on account of the high effectiveness of a magnetic field in 
exerting stresses over large distances, independent of the nature of the fluid in 
which is embedded (or even across vacuum). 

The 'memory-effect' of magnetic fields in HMD, i.e. their dependence on 
the history of fluid motions, allows large field strengths to build up over time, 
resulting in torques far larger than can be produced, for example, by the in- 
stantaneous overturning action of a hydrodynamic instability. But in principle 
this same effect also makes magnetic stresses dependent on initial conditions as- 
sumed for the magnetic field of a (differentially rotating) star. In Spruit (2004) 
I have argued that this dependence in fact leads to a relatively straightforward 
'bifurcation' of possibilities. 

If the initial field is strong (in some sense that can be quantified), inter- 
nal magnetic stresses eliminate differential rotation. Such strong fields have 
their own internal instabilities, causing the field to decay towards a minimum 
energy state, which depends on the degree to which magnetic helicity is con- 
served (Braithwaite and Spruit 2004, Braithwaite and Nordlund 2006), resulting 
in a stable uniformly rotating end state such as observed in magnetic A-stars. 
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Weak initial fields, on the other hand are wound up by the differential rotation 
present in the initial state (or developing as a consequence of stellar evolution), 
producing an azimuthal magnetic field (with respect to the rotation axis) which 
increases linearly with time until it becomes unstable under it own internal 
forces. In Spruit (2002, hereafter S02) I have outlined how such instability 
can feed a field amplification cycle producing a small-scale time dependent 'dy- 
namo'. The steady amplitude of this amplification process is then independent 
of the initial weak field assumed, and depends only on local conditions of differ- 
ential rotation, thermal and compositional stratification, magnetic and thermal 
diffusion. 

The estimates made in S02 depend critically on the conditions for instability 
of an azimuthal magnetic field, since it is this instability that closes the field- 
amplification loop. Fortunately, at least the linear stability conditions for an 
azimuthal magnetic field can be analyzed in great detail since the instabilities 
are intrinsically local in a meridional plane (though global in the azimuthal di- 
rection). Instability conditions were derived in an appendix in Spruit (1999, 
hereafter S99), on the basis of a general dispersion relation derived by Ache- 
son (1978). They include the effects of stabilizing thermal stratification and 
gradients in composition, thermal and magnetic diffusion. 

1.2 Observational tests 

In a star like the Sun the composition is still nearly homogeneous, with exception 
of the Hydrogen burning core and the effect of gravitational settling of Helium 
below the convection zone. In this case, the stabilizing effect of the thermal 
stratification, though strongly relaxed by thermal diffusion, is the dominant 
influence on instability. The slow and nearly uniform rotation of the radiative 
interior of the Sun is compatible with angular momentum transport by magnetic 
fields according to the dynamo estimate of S02. 

The internal rotation of the Sun is a key test case for any theory of angular 
momentum transport in stars, first of all because it rotates so slowly that most 
hydrodynamic transport mechanism fail to explain the observed uniformity of 
its internal rotation. This led to the early suggestion (Spruit et al. 1983) that 
magnetic fields are the dominant transport mechanism in the radiative interior 
of the Sun. A possible exception is transport by internal gravity waves excited 
by the convection zone (Charbonnel and Talon 2005). Secondly, the Sun is also 
the only star for which a value for the spindown torque, independent of the 
rotation rate, exists (from measurements in the solar wind). Some prescriptions 
for internal torques commonly used in stellar evolution codes ('with rotation') 
fail to reproduce the internal rotation of the Sun by large factors. 

Observations of rotation rates in various stages of stellar evolution give a 
confusing picture, however, and it seems unlikely that a single mechanism can 
fit all. Initial spins of pulsars around 10 ms are consistent with the magnetic 
torques of S02 (Heger et al. 2005). The collapsar model for gamma-ray bursts, 
on the other hand, requires much larger amounts of angular momentum to 
survive in a pre-supernova core. Thus any single, deterministic mechanism for 
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angular momentum transport can not fit both pulsar rotation and GRB. 

The rotation of white dwarfs can be explored as an indication of rotation 
rates in the cores of giants, and hence indirectly also of internal angular momen- 
tum transport during the giant phase. The observed rotation of WD scatters 
widely, however. In Spruit (1998) I have argued that the rotation of WD is in 
fact not a remnant of the initial angular momentum of the star, but a result of 
small random asymmetries in the mass loss on the AGB. A similar effect ('birth 
kicks') has been proposed for the rotation of neutron stars (Spruit and Phinncy 
1998). If it contributes significantly to pulsar rotation, the contribution from 
the initial angular momentum must be less, implying that the prescription in 
S02 in fact still underestimates the internal torques in a star. The presence 
of a kick contribution to pulsar rotation can actually be tested observationally 
by correlation of spin axes and proper motions of pulsars (Spruit and Phinncy 
1998, Wang et al. 2006). 

The rotation rates of horizontal branch stars are a potential source of in- 
formation on angular momentum transport in stars on the giant branch (de 
Medeiros 2004 and references therein). The observation that the more compact 
blue HB stars tend to rotate more slowly than red HB stars, however, is a chal- 
lenge for any mechanism for the angular momentum transport in their giant 
progenitors. 

Inferences about internal torques in stars have been made from patterns of 
rotation in young clusters of different ages. Differential rotation between core 
and envelope surviving for 10 yr fits these data (Dennisenkov and Pinson- 
ncault 2006 and references therein). The strong internal coupling implied by 
the dynamo estimate of S02, on the other hand, is not compatible with this 
interpretation. Dennisenkov and Pinsonneault (2006) propose a modification of 
the stability conditions in S99, resulting in lower magnetic torques. The mathe- 
matics of this modification contains an error, however. This is explained below 
in a paragraph following eq. (JSJ) . 

Calculations of the internal rotation of evolving stars including the magnetic 
torque prescriptions in S02 have been done by Heger et al. (2005), Maeder 
and Meynet (2005), and Yoon and Langer (2006). They show that regions of 
homogeneous composition quickly become uniformly rotating, with gradients in 
rotation rate accumulating in the composition gradients. Angular momentum 
transport across /^-gradients is thus the critical factor for the internal rotation 
of evolved stars. For this reason the instabilities in a composition gradient are 
scrutinized a bit closer here than was done in the appendix of S99. This is done 
in section [3J 

2 Instabilities without a //-gradient 

In this section the mathematics leading to the stability conditions in S99 is 
spelled out in a bit more detail, and expanded to include a discussion of the 
case when the magnetic frequency is larger than the rotation rate. 

Near the pole, and in the limit va <C c s Acheson's dispersion reduces to eq. 
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(A8) in S99, and can be written as 



w|(fim + w) + ^- " a j [wfw + iris 2 ) - ui] - 

2fim + s^(p + !)] I wm ( w + i? ? s2 ) + w i] = 0. (i) 

Here w is the (complex) frequency of the perturbation, I, n the horizontal and 
vertical components of the wave number, m the azimuthal order, s 2 = l 2 + n 2 , 
p the index of the variation of the field strength with distance w from the axis 
p = d In B/d In zu, and wa the characteristic magnetic frequency wa — v^/uj. N 
is the buoyancy frequency, 77 the magnetic diffusivity, k the thermal diffusivity. 
Viscosity has been neglected. 

In all of the following, I assume to > 0. Axisymmetric perturbations have 
to be analyzed separately from but they are unlikely to be very interesting. 

Equation Q is of fifth order in u, hence analytically tractable only in limiting 
cases. Fortunately, the parameter space most likely to be relevant in stars has 
an ordering of time scales: 

N > fl > uja- (2) 

[The latter of these inequalities is not immediately obvious; this is discussed 
separately in section l2~Tl below.] To see what the limit Q ^> wa implies, introduce 
the following dimensionless quantities: 

Q s 2 riQ s 2 kQ 
a = uj—^-, h= — 5— ; k= 5-, (3) 

(note difference of definition of h,k compared with S99). After multiplication 
by (a + ih){a + ik) eq. QJthen reads (without assuming e to be small yet): 

(to + ae )(p — l)(a + ih)(a + ik) 
+^-[ae 2 (a + ih)(a + ik) - (a + ik) - Q(a + ih)][ae 2 (a + ih) - 1] 

-[2m(a + ih)+p + l][m(a + ih) + l](a + ik) = 0. (5) 

All physical effects (buoyancy, magnetic energy, thermal and magnetic diffu- 
sion), contribute to the solution of this equation if all terms and factors, includ- 
ing the frequency a, are of the same order. One of the things this implies is 
that Q must be O(l), or 

l/n ~ N/uj A < 1, (6) 

that is to say, the radial scale of the mode must be much smaller than its 
horizontal scale. If l/n is not small, the term involving Q dominates, and must 
be balanced in the first square bracket in JSJ by a suitably large u>. The solutions 
for such wavenumbers are thus buoyancy waves (gravity waves), uj 2 — (l 2 /s 2 )N 2 . 
For somewhat smaller values of Q, but still 3> I, a balance is possible with 
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frequencies around the rotation rate Q; these rotationally modified are buoyancy 
waves are called inertial-gravity waves. With the ordering of time scales © they 
are unaffected by the magnetic field. To find the modes in which the magnetic 
field plays a role, one must have Q ~ 0(f). This is the formal justification of 
the heuristic derivation of the minimum vertical wavenumber in S99 (eq. 44). 
Note that it does not involve the rotation rate tt. 

Magnetic and thermal diffusion become important for wavenumbers such 
that h,k ~ 0(1), respectively, or 

. ,2 , ,2 

The typical value of a for a solution of J^J is of order unity, and thus the growth 
rate a = — Im(u;) of an instability of order ~ lu\/Q. Hence we can interpret (J2J 
as stating that the effect of magnetic diffusion on instability becomes important 
when the growth time is of the same order as the diffusive time scale l/(n 2 r]) 
(and similar for thermal diffusion k). The rotation rate therefore does have 
an effect on the high end of the wavenumber range, where diffusion becomes 
important. This agrees with the heuristic argument in S02. 

Now seeking solutions for the case c < 1, but keeping all other quantities 
in J5J of order unity we can expand © in powers of e 2 . It appears only in the 
combination ae 2 , with the result that at the lowest order in e, JSJ is only of 
third order in a (eq. A18 in S99): 

m(p - l)(a + ih)(a + ik) + ^[a(l + Q) + ih + ik] 

-[2m(a + ih)+p+ l][m((a + ih) + l](a + ik) = 0. (8) 

The two remaining roots of the original fifth order equation have disappeared 
because they do not scale as u> ~ uj^/Q. 

The procedure to derive instability conditions from this is now as follows. 
Since the coefficients in (JHJ are complex, the solutions a are complex as well. 
Roots with Im(a) > decay as exp[— Im(a)] exp[iRe(a)], hence correspond 
to stability. Roots with Im(a) < correspond to instability. The boundary 
between stable and unstable modes is at Im(a) = 0. To find the instability 
condition, we thus set lm(a) = and take Re and Im parts of (JSJ), yielding 
two equations for the real quantity a that have to be satisfied simultaneously. 
Together, they determine not only the value of a, but also a relation between the 
parameters that has to be satisfied; this relation defines the boundary between 
stable and unstable roots. 

In their discussion of the dispersion relation (J5J Dennisenkov and Pinson- 
neault (2006) misinterpret this procedure. By assuming a to be real (as above) 
and eliminating the highest power of a between the real and imaginary parts 
of Q they obtain an equation of lower order. They proceed to obtain complex 
solutions of this equation, apparently motivated by the fact that complex values 
of the mode frequency indicate instability. This is incorrect since the separation 
into real and imaginary parts used holds only for Im(a = 0), and it ignores the 
fact that there are two equations to be satisfied. 
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The imaginary part of © can be written as: 



- 2m 2 a 2 - 4ma - 4^(m 2 a 2 + ma) + ^-(Qt + 1) - (p + 1) + 2m 2 h 2 = 0. (9) 
k 2 k 

It is not necessary to consider the general case where both diffusivities are 
similar. In a stellar interior the thermal diffusivity is generally much larger 
than the magnetic diffusivity. In taking the limit h/k — 777/K <C 1 in J^J), the 
term involving Q must be kept, however. It represents the stabilizing effect 
of the stratification, proportional to the square of the aspect ratio l/s of the 
perturbations. The condition Qi]/k ~ 0(1) determines the range of aspect ratios 
that is relevant for the effect of the stable stratification on the instability. Due 
to the factor 77/ ' k multiplying Q compared with the adiabatic case 77 = k = 0, we 
can already conclude that the range of unstable aspect ratios is now larger than 
in the adiabatic case. Effectively, thermal diffusion reduces the buoyancy iV 2 by 
a factor tj/k. This is very similar to the effect it has on purely hydrodynamic 
instabilities (Zahn 1974). Thus © becomes: 

-2m 2 a 2 ~4ma+^-(Q^ + l)-(p+l) + 2m 2 h 2 = 0. (n/n < 1) (10) 
2 k 

This can be written as 

2 

2 - 2(ma + l) 2 + ^- - (p + 1) + f + g 2 = 0, (11) 



where 



„9 m 2 n m 2 l 2 N 2 n 9 9 /x 

/ = -^Q^z = -iriCTr* = 2 ™ & ■ ( 12 ) 



2 K 

The real part of (JSJ) yields: 

2 

- 2m 2 a 3 - 4m 2 o: 2 + a[2m 2 h(h + 2k) - (p + 1) + ^(Q + 1)] + Amhk = 0. (13) 

It is not immediately clear from this equation how the limit h/k <C 1 is to be 
taken, but multiplying (|10J) by a, subtracting this from jjE} , and multiplying 
the result by /i/fc yields: 

a[ g-(l--)+4m 2 / l 2 ]+4m/t 2 ^0. (»//«< 1) (14) 

2 k k 

Since Qh/k ~ O(l) the limit /i/fc — > thus yields 

m Q = -2j 2 /(/ 2 + S 2 ). (15) 

This implies that 

2-2(rna + l) 2 > 0. (16) 

Hence all terms in 111 U are positive definite except the term involving p. A 
necessary condition for instability is thus that a = p + 1 — m 2 /2 > 0. This is 
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most easily satisfied for m = 1 (recall that the case m — has to be treated 
separately). For azimuthal fields generated by winding up near the pole, p = 1, 
hence a = 3/2 and the condition on p is satisfied for m = 1. The other necessary 
conditions are: 

2 - 2(ma + l) 2 < a and f 2 < a and g 2 < a. (17) 

Sufficient, on the other hand is the set of conditions 

2-2(ma + l) 2 < a/3 and / 2 < a/3 and g 2 < a/3. (18) 

The necessary are close enough to the sufficient conditions. The gap between 
the two can be narrowed with a more detailed analysis but the present approx- 
imation is sufficient for application to the dynamo estimates in S02. One then 
finds that, up to a factor of order unity, a necessary and sufficient condition for 
instability is (S99 eq. A29:) 

I I S 

lo%> l -^prl 2 N 2 ^Vt, (k>77, ^>cj a , p>-l/2, m>0). (19) 
a"' k 

One may wonder what has happened to the two roots that have disappeared 
from the original 5th order equation, and whether those could become unstable 
as well. In the limit 3> loa, these roots turn out to scale as ui ~ fl. Under 
this assumption, and expanding in e as above, the dispersion reduces to 

72 

lu 2 - —N 2 - AQ 2 = 0. (20) 

The magnetic field has thus disappeared from the equation, showing that the 
two remaining roots purely hydrodynamic. They describe inertial-gravity waves, 
and are unstable only for N 2 < 0, i.e. in convective zones. 

2.1 The 'slowly rotating' case 

One may wonder if the assumption f2 3> wa is always justified, and how the 
instability would be different in the opposite case. This question is also moti- 
vated by the fact that the dynamo argument in S02 does not lead to a properly 
defined amplitude for the magnetic field generated in this case. 

It is now more convenient to write the dispersion relation in terms of a 
dimensionless growth rate a (which, for the moment, is still a complex number): 

lu = —iauJA- (21) 

Setting 

5 = Q/uj a , H = t]s 2 /lo a , K = ks 2 /-/lj a , (22) 
the dispersion relation Q can be written as 

( ff + imS) - + _i_ + ^}[a(a - H) + 1] 

+(2imS- ^jj)[im5(a - H) - 1] = 0. (23) 
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To lowest order in S, the coefficients of this equation are real. The condition 
for exponential instability is found from the sign of the constant term after 
multiplying by (a — H)(a — K), and is: 

P>-l + ^(l + Q-). (24) 

Since Q can be made arbitrarily small by choice of the wavenumbers, there are 
unstable modes when p > — 1 + m 2 /2. This condition is as before in the rapidly 
rotating case: the (logarithmic) gradient in field strength required for instability 
is the same. At this lowest order in D,/u>a, there is no threshold for instability 
for an azimthal magnetic field in a nonrotating star; the only requirement is 
that the gradient index p be large enough. To find the minimum field strength 
for instability, the next order in the small quantity 6 has to be considered. With 
a bit of analysis one finds 

m 

p> -l + — +8mS, (25) 

i.e. a small correction for 6 <C 1. More significantly one finds that there is now 
a minimum field strength for instability: 

7r*(£>"^>" 4 <^v>^ («»"> <*» 

This is the same condition as for the opposite limiting case wa/^ *C 1 (see eq. 
(A29) in S99, noting that the equivalent eq. (49) in S99 contains a typo, as 
explained in a footnote in S02). Thus (|26|l happens to apply irrespective of fi. 



2.2 Relevance of the 'slowly rotating' case 

Next, we can ask when the case fi <C ua is likely to apply. For this, we need 
an estimate of the amplitude of the magnetic field. If it is generated by the 
dynamo process envisaged in S02, and if the estimate for its amplitude given 
there is valid, the field strength in the rapidly rotating case would be given by 

(S02 eqs. 17, 19), where q is the dimensionless rotation gradient q = din fi/d In 7-. 
This estimate applies if the buoyancy frequency N is dominantly due to the ther- 
mal stratification (in a sense that can be defined more precisely). If, instead, 
the buoyancy due to a gradient of composition is more important ('case 0' in 
S02), the field strength is given by 

Q-^AT (28) 

where N = is now the buoyancy frequency due to the stratification of 
composition. For the slowly rotating case uj^/SI > 1 to apply, the field strength 
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generated by the dynamo has to be high enough. From eqs. (|2"7|) . l(2*%)l it follows 
that this is possible only for steep rotation gradients g> 1. Such steep gradients 
can occur in two cases. 

In normal (slow) stellar evolution, the high effectiveness of the magnetic 
torques keeps the level of differential rotation low wherever the composition is 
uniform. The gradients in rotation rate then concentrate in the /^-gradients (e.g. 
Heger et al. 2005). Thus, the most likely place where the slowly rotating case 
might apply is when the composition gradient determines the buoyancy. Eq. 
<(2*5|) then shows that for id a to be larger than f2, one must have q > N/Q, or 

rdfl/dr > N. (29) 

Up to a numerical factor of order unity, this is the Richardson condition for 
(adiabatic) hydrodynamic instability in a stable stratification. The 'slowly ro- 
tating' case n <C wa is thus of limited relevance, since where it applies, ordinary 
hydrodynamic instabilities will already be effective. 

In the above I have assumed that the evolution of the rotation profile in 
the star is slow enough that the dynamo process is in a quasistationary state, 
so that winding-up of the field by differential rotation is balanced by its dacay 
due to magnetic instabilities. If, instead, the rotation profile evolves faster than 
the time scale on which the dynamo attains its equilibrium amplitude, the field 
can increase to larger values, solely by the rapid winding-up of field lines. This 
might occur in pre-supernova cores. 

3 Instabilities in a /i-gradient 

If the buoyancy is due to a composition gradient, it can not be reduced by 
thermal diffusion. The only relevant diffusion affecting the buoyancy is then the 
diffusivity, A, of the ions responsible for the stability of the stratification. This 
diffusivity is generally several orders of magnitude smaller than the magnetic 
diffusivity (which is mediated by the electrons). This case, A <C i], is thus 
equivalent to setting K/7 = A <C rj in the (Q, while interpreting N as the 
buoyancy frequency iV M = (gd r \n [i) 1 / 2 due to the ^-gradient alone. 

In this case, the agent driving the instability (the magnetic configuration) 
has the higher diffusivity, the composition gradient causing the stabilizing strat- 
ification the lower diffusivity of the two. As in the case of other double-diffusive 
instabilities (e.g. Turner 1980), we should expect the existence of an intrinsi- 
cally oscillatory instability (i.e. even when rotation is absent). An example 
is the case of hot salty water under cooler fresh water, such that the overall 
density stratification is stable. Internal gravity waves in this stratification can 
grow on wavelengths such that the thermal diffusion time is comparable with 
the wave frequency. The astrophysical equivalent is 'semiconvection', where the 
corresponding linearly unstable modes are called Kato oscillations (Kato 1966). 
Such modes should also be expected in the present case. 

The analysis is a bit more complicated than the case k ^> r\. To illustrate 
the nature and stability conditions for such modes, consider the simplified case 
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n = 0. As in the case of Kato oscillations, we expect the unstable oscillations 
to have frequencies near the buoyancy frequency F of a wave with vertical and 
horizontal wavenumber components n, I. Introduce the scalings and notation 

/ wi to ris 2 As 2 

F = -N, -± = f, -=,, h=*-, - = k (30) 

(noting that the symbols h, k have been recycled from previous usage above). 
The dispersion relation Q then becomes 

m 2 f 2 1 f 4 

The conditions for marginal overstability are found by assuming v real, and 
taking the Re,Im parts of this equation. The real part has a common factor v 
and after division by this factor reads 

f 2 (p l)(v 2 hk) + ^-hk-f- l]{v 2 f) 

9 

-^-\v 2 {h + k)-fk-h\h-f{y+\) =0. (32) 
This is a quadratic in v 2 . The Im part yields 

/V( P - i)(i + \) + ^(1 + k j)-f\- 1](- 2 - 1 2 ) 

-^ 2 -hk-f 2 -iy-f( P +l)^ =0, (33) 

also a quadratic in v 2 . 

So far, no assumptions about large and small have been made. The quantity 
/ measures the strength of the magnetic field relative to the stratification. For 
/ ^ 1, the field is strong enough to overcome the stable stratification irrespective 
of diffusive effects. Since we are interested in just these effects, we may assume 
the limiting case / C 1. By assumption, we also have A/77 = k/h <S 1. Not 
making further assumptions about the relative smallncss of / and k/h, and 
choosing wavenumbers such that h 2 is also small, we set 

f = e, k/h ~ 0(e), h 2 ~0(e). (34) 

To order zero in e, both loL'l) and (|33|l then yield v 2 = 1, so 

v 2 = l + 0{e). (35) 

The mode frequency is thus known now, and confirms expectation that it is the 
frequency of a buoyancy wave. We are thus looking for the conditions under 
which a buoyancy wave can become overstable. To find the stability condition, 
the next order in e has to be considered. It is found by eliminating the term in 
v A between (|3^|) and Expanding the result to lowest order in e then yields 

m 2 k , , . 

^-l--(l + ^) = 0, (36) 
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as the condition that holds at marginal oscillatory instability. Evidently, a 
necessary condition for instability is thus 



P>1 + — . (37) 

Setting a = p — 1 — m 2 /2, 1|36|1 becomes, in dimensional quantities: 

o I 2 m 2 A . A , . 

u 2 A = —N 2 — - (-«1, o>0). (38) 
n z 2a r) rj 

Though we formally do not know from this analysis on which side of equality 
the instability lives, it clear that l|38[) actually defines a minimum field strength 
for instability, since the field is the ingredient driving the instability. 

Condition (|38() is similar to the opposite case n/n <C 1, but with r\ and n 
reversed. This is characteristic of double diffusive instabilities in general. Note, 
however, that the requirement on the field gradient index p is different from the 
case tj/k <C 1, where we had p > — 1 + m 2 /2 instead. The oscillatory instability 
in a fi- gradient thus requires a significantly steeper field gradient. In particular 
a field would up by differential rotation near the pole, p = 1, is stable. Also 
note that the terms h 2 / f 2 have canceled in the result l|38|l . This has the effect 
that there is actually no minimum wavenumber for instability, and hence no 
minimum field strength cither, since l/n can be taken arbitrarily small in (|38|) . 
This differs from the case r\j k, <C 1, where at the same order there is a minimum 
field strength (cf eq. I26fl . It must be noted, however, that the absence of a 
minimum field strength is probably an artefact of the inviscid approximation 
used here. Viscosity is likely to define a minimum field strength. In view of the 
uncertain relevance of this instability, and the likelihood that linear stability 
gives a misleading impression anyway in this case (see section fc.lfl . this is not 
explored further here. 



3.1 layering 

If a /i-gradient is initially strong enough to prevent instability, evolution of the 
star will tend to increase the gradient in rotation rate across it until instability 
happens, either by the process described above, or the simpler form of instability 
on which the dynamo estimates of S02 are based. Once this happens, however, 
it is quite possible that the subsequent nonlinear development will lead to layer 
formation: the gradient in [i breaks up into a series of steps, each with (nearly) 
uniform composition and rotation rate, separated by steep gradients in [x and 
Q. If this happens, the angular momentum transport can potentially increase 
strongly. 

To see this, imagine replacing the smooth gradient in O and /i by a series of 
steps of thickness d, with boundary layers of some width between the steps. In- 
side a step, the stabilization by a /^-gradient has disappeared, so that instability 
and dynamo action can take place as if there were no gradient in composition at 
all. If the instability time scale is er -1 = fi/^A (see section [21 above) , there will 
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be magnetic boundary layers of width 5 = {rj/a) 1 / 2 bounding each step. Across 
these very steep boundary layers, the exchange of magnetic flux is entirely by 
magnetic diffusion rather than fluid motion. 

To set up such a layered structure, energy has to be spent to overcome the 
stable composition gradient, but once it is in place, it can persist. The energy 
required to overturn the stratification into layers of thickness d is, per unit 
mass, e w d 2 N 2 (equal to the kinetic energy of a vertical buoyancy oscillation 
of amplitude a). If the composition gradient extends over a radial distance 
D, this energy is smaller by a factor (d/D) 2 than the energy per unit mass 
needed to overturn the stratification globally. This small amount of energy 
might be supplied by, for example, gravity waves from a nearby convective zone, 
or else by magnetic instability when the gradient in rotation rate has become 
strong enough in the course of evolution of the star. Once layering sets in, the 
stabilizing effect of the composition gradient is reduced. This is quite analogous 
to the case of double diffusive instabilities in geophysics and in the laboratory 
(c.f. Turner 1980, 1985). 

An astrophysical equivalent of this nonlinear, subcritical behavior of insta- 
bilities in a composition gradient is semiconvection. Once a layer within a 
composition gradient has been overturned, be it by a numerical artefact or by 
an actual physical process, the mixed state of the layer will be maintained by 
convection. This is often observed in semiconvective regions of stellar models, 
even if layer formation (Spruit 1992) is not taken into account explicitly. 

The same may be taking place in stellar evolution calculations which include 
magnetic torques according to S02. In the results of Heger et al. (2005), for ex- 
ample, there is evidence of spiky structure in the /i-gradients on length scales 
comparable to the numerical resolution. To assess its consequences quantita- 
tively, it would be important to find ways of taking the layering process into 
account in a more controlled way. 
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